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We discuss the issue of measuring the mean position (center-of-mass) of a group of bosonic or 
fermionic quantum particles, including particle number fluctuations. We introduce a standard quan- 
tum limit for these measurements at ultra-low temperatures, and discuss this limit in the context 
of both photons and ultra-cold atoms. In the case of fermions, we present evidence that the Pauli 
exclusion principle has a strongly beneficial effect, giving rise to a 1/A'' scaling in the position 
standard-deviation - as opposed to a scaling for bosons. The difference between the ac- 

tual mean-position fluctuation and this limit is evidence for quantum wave-packet spreading in the 
center-of-mass. This macroscopic quantum effect cannot be readily observed for non-interacting 
particles, due to classical pulse broadening. For this reason, we also study the evolution of photonic 
and matter-wave solitons, where classical dispersion is suppressed. In the photonic case, we show 
that the intrinsic quantum diffusion of the mean position can contribute significantly to uncertainties 
in soliton pulse arrival times. We also discuss ways in which the relatively long lifetimes of attractive 
bosons in matter-wave solitons may be used to demonstrate quantum interference between massive 
objects composed of thousands of particles. 

PACS numbers: 03.75.Pp, 42.50.Dv 



I. INTRODUCTION 

The topic of mesoscopic quantum effects is of great 
current interest, especially in photonic or in ultra-cold 
atomic systems where a large decoupling from the envi- 
ronment is possible. An important degree of freedom is 
the average or center-of-mass position of a large num- 
ber of particles. It is now possible to observe quantum 
diffraction of the center-of-mass in molecules like Cgo and 
related fuUerenes jlj] . Larger physical systems would pro- 
vide an even a stronger test of quantum theoretical pre- 
dictions. As well as testing quantum theory in new, meso- 
scopic regimes, these types of experiment have potential 
applications to novel sensors. For example, ultra-precise 
measurements of position may be useful in measuring 
gravitational interactions, or other forces that couple to 
conserved quantities. 

In this paper, we wish to examine the quantum limits 
to center-of-mass position fluctuations. We further dis- 
cuss these limits in the context of mesoscopic systems of 
interacting particles, atoms or molecules. This is an im- 
portant issue in quantum-optical or ultra-cold atom en- 
vironments, where center-of-mass motion places a limit 
on coherence properties of photonic or atom lasers. A 
localized system of particles features a dual particle and 
wave nature, that is, it has conjugate observables of mo- 
mentum and position, as well as conjugate observables 
of number and phase. As pointed out by Landau and 
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Peierls[2], free particle momentum is a QND observable. 
In fact, the momentum of a quantum soliton can be non- 
destructively measured via the position of a probe soliton 
after their collision^. The quantum position fluctuation 
increases with propagation distance due to an initial mo- 
mentum uncertainty. This places a fundamental limit on 
momentum QND measurements. 

We generalize the concept of the standard quantum 
limit of a mean position measurement ^, ^ to any kind 
of quantum field. The standard quantum limit is defined 
here as the position uncertainty of a many-body ground- 
state for non-interacting particles with a given density 
distribution and total momentum. For massive particles, 
this corresponds to the Heisenberg-limited uncertainty of 
the center-of-mass of a group of non-interacting particles 
in an external potential, at zero temperature. Intrigu- 
ingly, we find completely different scaling laws for dif- 
ferent particle statistics: while bosons have a position 
uncertainty (standard deviation) that scales as 1/^/N, 
fermionic center-of-mass uncertainties scale as 1/N for 
N particles. This is caused by the particle correlations 
induced by Fermi statistics. The scaling indicates that 
ultra-cold fermions are the preferred system for ultra- 
sensitive measurements of center-of mass. 

A necessary requirement for quantum- limited measure- 
ments is an extremely low level of fluctuations and noise, 
which makes laser pulses and ultra-cold atoms the most 
reasonable choices at present. As an example, the lo- 
calized solitons of the attractive one-dimensional Bose 
gas can now be observed both in quantum-optical and 
in atom-optical environments. As well as technologi- 
cal applications in communications, these types of soli- 
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ton show intrinsic quantum effects. That is, the ef- 
fects of quantum phase diffusion (squeezing) [1, 0, 0] have 
been observed in experiments [1, M, E^- Quantum cor- 
relations estabhshed by soUton coUision (quantum non- 
demohtion measurement) [3] have also been demonstrated 
in experiment pH . Even though they involve up to 10^ 
particles, these effects simply do not occur in classical 
soliton theory. 

Phase squeezing or phase QND measurement does not 
lead to center-of-mass changes, which are especially in- 
teresting for massive particles, as this degree of freedom 
is coupled directly to the gravity field. Due to prob- 
lems in quantizing gravity, there have been suggestions 
that gravitational effects may lead to wave-packet col- 
lapse and/or dissipation [H, [ll, [H, [H, [IJ- While this 
remains speculative, it is clearly an area of quantum me- 
chanics where there are no existing tests. To exclude such 
theories, one would need to violate a classical inequality 
involving mesoscopic numbers of massive particles, whose 
center-of-mass was in some type of quantum superposi- 
tion statefl^. In real experiments, weak localization 
can take place due to interactions with the environment. 
It is at least interesting to investigate how this may occur, 
as a first step towards testing more fundamental types of 
localization. We note that steps in this direction have 
been taken recently for non-interacting, massless parti- 
cles, via position observations with photons, at the stan- 
dard quantum limit or better [20l. [2]|. 

Quantum wave-packet spreading of the center-of-mass 
is an intrinsic quantum prediction for untrapped parti- 
cles at all particle numbers. However, this is masked 
by the single-particle effects analogous to classical pulse 
broadening for an optical pulse in a linear dispersive 
medium. A soliton or quantum bound-state can suppress 
this single-particle or classical pulse broadening, and al- 
lows the intrinsic center-of-mass quantum spreading to 
be observed [[^, in a similar way to the known quantum 
noise effects found in amplified solitons|23,]. We find that 
this is an exactly soluble problem for any initial quantum 
state. No linearization or decorrelation approximation 
is needed, which allows us to examine quantum wave- 
packet spreading over arbitrary distance scales. In or- 
der to distinguish classical pulse broadening from quan- 
tum wave-packet spreading, we calculate the standard 
quantum limit of this measurement for an initially coher- 
ent soliton pulse. This requires knowledge of the time- 
evolution of the density distribution for an interacting 
many-body system, which is nontrivial, and requires nu- 
merical solutions. 



II. CENTER OF MASS OPERATORS 

We start by noting that while a centroid exists for any 
localized, measurable quantity, it is most interesting for 
a conserved quantity. For any physical system with a 
conserved current four-density j = (joj), it is possible to 
define a center-of-charge (or mass, energy, etc) relative 



to the conserved quantities 

/ j^(x)d^x 



(1) 

In a quantum state that is an eigenstate of the conserved 
charge with Jq > 0, defined on a D- dimensional space, 
the definition of the centroid is: 

After partial integration, the conservation law. 



dt 



(3) 



then implies that the centroid position translates with a 
velocity given by the spatial part of the conserved cur- 
rent: 



J 

Jo 



(4) 



While the concept of measuring the center of mass 
position is well-known, there are a number of issues in- 
volved. The most obvious is that there is a difference 
between the center-of-mass and the average position of 
a group of non-identical particles. This difference van- 
ishes for non-relativistic particles all of the same mass: 
for example, isotopically pure, ultra-cold atoms. How- 
ever, there is a real difference in the case of particles of 
different masses, and for extremely relativistic particles 
and photons, where the effective mass depends on the 
momentum. Here the average position differs from the 
center-of-mass. Another subtlety which this paper treats 
in some detail, occurs when the system is in a quantum 
mixture of different particle numbers: the effects of this 
are treated in the remainder of this section. 



A. Massive fields 

We first consider non-relativistic massive quantum 
fields '^i (x) , where each component has an equal 
mass. The fundamental particle statistics can be either 
fermionic or bosonic. As usual for quantum fields, the 
commutators are: 



*,(t,x),^'Ut,x') =5y53(^„^/) 



Introducing the particle density 



(5) 



n(x)=^n,(x)=5]vl;t(x)^,(x), (6) 



we can define the total particle number 

N ^ ( d^xh (x) , (7) 
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which is exactly conserved in the absence of dissipative 
processes hke absorption or amplification. In this paper, 
we will consider cases where the system has an uncer- 
tainty in N, but we will exclude states with zero parti- 
cle number. In terms of operational measurements with 
small particle numbers, this may require post-selection 
to exclude states of this type - since of course the center- 
of-mass is undefined in such cases. 

In order to define center-of-mass or mean position for 
systems with particle number fluctuations, several defi- 
nitions are possible, depending on measurement proce- 
dures. There is more than one measurement procedure 
possible, since the quantum state need not be an eigen- 
state of the number operator. Hence, different number 
states could have an arbitrary relative weighting. 

In the absence of external potentials, systems of inter- 
acting particles also have an invariant quantity due to 
translational invariance, which is the total momentum 
(in D dimensions): 



h 



^ (x) V*, (x) - (x) (x 



(8) 



Here we use the Einstein summation convention for re- 
peated indices. We use a capital letter to emphasize that 
this is an extensive quantity, proportional to the number 
of particles. 



that: 



£ = E + iB 

= V xA 
= £+ + £^ . 

It is easily checked from the complex form of Maxwell's 
equations that id£a/dt = V x f^and idAa/dt = V x Aa- 
The photon density is then: 

n(x) = ^n.(x) = ^:(£t.X+4-it) (9) 



with a corresponding conserved current of 



J = 



xAa + Aix£a 



(10) 



We note that with these definitions, the photon den- 
sity is not positive definite in small volumes and time- 
intervals. However, the total photon number N is well- 
defined and positive definite, so that a position centroid is 
well-defined for this conserved quantity. For the remain- 
der of this paper we will focus on non-relativistic massive 
particles for definiteness, while noting that many results 
that only depend on the existence of a conserved current 
will hold in the general case of a photon field. 



Massless fields 



C. Intensive position 



Massless fields - photons - are the most commonly used 
particles where measurements are made at the quantum 
level. In this case, we can distinguish two different im- 
portant cases. 

• Narrow band fields in one-dimensional wave-guides 
or fibers experience dispersion, which causes them 
to behave as massive particles. In a similar way, the 
paraxial approximation for beams allows diffractive 
behavior to be treated in an effective mass approx- 
imation. These cases can be handled in the same 
way as the treatment given above. 

• More fundamental issues arise in free-space, where 
there is a long-standing fundamental problem of 
how to define photon position. In addition, as there 
is no photon mass, even the idea of center-of-mass 
requires care. Instead, one must either use the con- 
cept of a center-of-energy, or of a center of photon- 
number. We shall focus on the latter concept here. 

To treat the free-space average position, we can use a 
conserved local density obtained from the dual symmetry 
properties of Maxwell equations [2^1 . We introduce £„, 
where a — ±1, which are the helicity components of the 
complex Maxwell fields in units where c = eg = iiq — 1, 
and the corresponding complex vector potentials Aa , so 



For quantum states that are eigenstates of number 
with > 0, the obvious definition is: 



■AN) _ 



J x.h (x) d^x 



N 



(11) 



This is undefined for the vacuum state, but is well-defined 
for all other number eigenstates. For the non-relativistic 
massive particle case, this obeys the expected commuta- 
tion relation with the total momentum operator: 



The presence of the vacuum in states composed of mix- 
tures or superpositions of number eigenstates introduces 
the requirement of deciding which relative weights to use 
for different particle numbers. The important fact to 
note is that the mean position of any vacuum compo- 
nent is undefined and, as such, its presence does not 
contribute to one's knowledge of the mean position. Ac- 
cordingly, as mentioned previously, we choose to work in 
a restricted Hilbert space which does not include the vac- 
uum state. Practically, this is equivalent to performing 
heralded particle number measurements where any null 
results are discarded. From a mathematical perspective, 
we work with a projected density operator p' = VpV, 
where V = (/ — |0)(0|). In order to extract position 
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information from linear superpositions or mixtures of 
number states, we introduce the center-of-mass position 
operator [iOl which is well defined in our projected Hilbert 
space: 



/ xn (x) d^x 



N 



(12) 



This intensive quantity also obeys the expected commu- 
tation relation with P: 



= ihSi 



(13) 



Thus, provided any measurement of the vacuum is dis- 
carded, X and P form a pair of canonical conjugate vari- 
ables with the dimensionless uncertainty relation 

{Ax'^){AP^)>^ (14) 
for states projected into our restricted Hilbert space. 



D. Extensive position 

Although X has the usual definition of the mean po- 
sition, it's operational denominator makes it is difficult 
to measure using many standard techniques, or even to 
represent in a normally ordered form. We can also in- 
troduce an extensive position operator, which produces 
a result proportional to the number of particles as well 
as their position. We define: 



xiV, 



(15) 



which is well-defined without projection in the Hilbert 
space. However, X and P do not form a pair of canonical 
conjugate variables, since 



(16) 



Through the uncertainty relation, we see that X is sen- 
sitive not only to positional fluctuations but also fluctu- 
ations in the total particle number or beam intensity: 



(17) 



Furthermore, as an extensive operator, X transforms in 
the following way. In a new reference frame S' where 
x' = X 4- Ax, we find that 



X' = X + AxTV 



(18) 



rather than just being translated by a c-number. 

We note that the extensive position does have a canon- 
ical conjugate partner, in terms of the intensive or mean 
momentum operator: 



. P 

with a corresponding commutator of: 



(19) 



(20) 



E. Quasi-intensive position operator 

We will finally define a quasi-intensive position opera- 
tor, which corresponds to the typical direct measurement 
procedures in optics. With this operator, the position is 
first measured with techniques that are extensive - but 
the final result is normalized by the average particle num- 
ber. This has the property that in an ensemble with vari- 
able numbers of particles, more weight is given to those 
measurements that involve the most particles. Here: 



(TV) 



(21) 



which can still be considered a useful measure of the mean 
position of a particle distribution, especially in the limit 
of large (N) where total number fluctuations are sup- 
pressed. Indeed, this is the variable that corresponds 
most closely to some earlier proposed operational mea- 
sures of pulse position, measured using local oscillator or 
homodyne techniques. It is interesting to note that the 
conjugate variable, the mean momentum, is the same 
as for an extensive position. This is analogous to the 
mean frequency of a laser pulse in a dispersive medium, 
which also corresponds closely to standard operational 
measurements used in laser physics. 

In order to measure fluctuations about (X) or (X), 
however, care must be taken due to the fact that X trans- 
forms in the following way under translation: 



X' = X + Aj 



(22) 



By considering fluctuations about x' = in a new refer- 
ence frame where x' = x — (X) , we define 



AX = X- (X)^ 

(N) 



(23) 



and thus minimize the number fluctuation contribution. 
Following from Eg. 1171 this obeys the uncertainty relation 



(AX?)(AP, 



> 



as the total number and momentum operators commute. 

In order to avoid complications arising from this issue, 
we will only consider isolated systems where it is always 
possible to flnd an inertial reference frame in which (X) — 
(P) = 0. 



III. QUANTUM LIMITS ON POSITION 
UNCERTAINTIES 



Intrinsic quantum uncertainty in the results of projec- 
tive measurements is perhaps the most striking difference 
between quantum mechanical observables and classical 
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variables. Here we seek to characterize this uncertainty 
in the center-of-mass (COM) measurement by introduc- 
ing a standard quantum limit to the variance in the dis- 
tribution of results. This is not a lower bound - there is 
none - but rather a natural limit that one can expect to 
achieve using standard cooling and/or stabilization meth- 
ods. The question of exactly which state to use to cal- 
culate such a characteristic COM is difficult to answer 
uniformly, especially as the states accessible to bosonic 
systems are different to those accessible to fermionic sys- 
tems. 

For a system with known density distribution and ex- 
change statistics, we define the standard quantum limit 
to the COM uncertainty to be the variance remaining 
when a non-interacting system of the same particles is 
reduced to zero temperature in an external potential that 
reproduces the given density distribution. It is important 
to note that this definition makes no restriction on the 
statistics of the characteristic system. However, we will 
find that the statistics do have a large effect on the way 
that the quantum limit scales for a given density. 

In the following analysis, we consider the variance in 
the intensive mean position operator, in the form: 

I Aip = ^ 1 d'^xd^y {Ax • Ayn(x)n(y)} . (24) 

Re-arranging this expression by using normal-ordering 
and commutators gives the result that: 



|AxP = J- (7V(t2+ : |AX|2 



(25) 



where we define the wavepacket variance operator 

^'-jI d^x{|Ax|2n(x)} (26) 

and the normally ordered COM variance 

: |AX| : = J d^xd^yAx- Ay 

x$]*l(y)*](x)^P,(x)^,(y). (27) 



Note that the exact operator ordering in the last expres- 
sion is important if the decomposition is to be correct for 
both Fermi and Bose field operators. 



A. Bosonic fields 

The simplest configuration for a system of degenerate 
bosons has all particles described by a single normalized 
mode function x(x). Thus, the characteristic states con- 
sidered here are formed using functions of the creation 

operator a^, = J d^x |x (x) 'I'^ (x)|. 



1. Number states 

Zero temperature bosonic number states can expressed 
in the following way: 



|iV) = -^(at)^|0) 



(28) 



Evaluating the COM variance of this state using the 
above decomposition is straightforward as the expecta- 
tion value of the normally ordered term : |Axp : van- 
ishes, leaving 



|Axp> 



N 



(29) 



where ct^ = (ct^) 



2. Cotierent states 



For bosons, the first term also vanishes for any co- 
herent state of the field operators, given our co- moving 
reference frame. Here we define a coherent state using 
a projection, as explained earlier, which projects out the 
vacuum state in order to correspond to a heralded or 
post-selected measurement. Such states can be formally 
defined as: 



-1/2 



E 

n=l 



|0), 



(30) 



This state can be viewed as a suitable reference state 
for quantum noise, in which the uncertainty in a position 
measurement is governed entirely by the spread in the 
particle distribution, |x(x)|^. After normalizing by the 
particle number, we find 



lAxI' 



(31) 



This is different from the following result, obtained by 
considering the variance in the quasi mean position op- 
erator X which for the coherent state above is given by 



(|AXp 



(32) 



From the graphical comparison shown as Fig. [T] it is 
clear that the standard quantum noise limit for the ex- 
tensive (quasi mean) position measurement is lower than 
that for the intensive position measurement. This can 
be understood as a consequence of the fact that an ex- 
tensive position measurement for a coherent state effec- 
tively weights all individual particle position measure- 
ments equally. This is advantageous for a coherent state 
- as in a laser pulse - where all particles arrive indepen- 
dently, and carry the same information. However, not all 



6 




0.0 



(|AXp)/a2 
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Figure 1: Deviation of the variance (relative to the wavepacket 
variance, a^) of the quasi intensive COM coordinate (dashed) 
from the true intensive COM coordinate (solid) for a heralded 
coherent state. 



pulses have the same number of independent particles. 
An intensive measurement weights all collective parti- 
cle position measurements equally, regardless of particle 
number. This gives less weight per particle measurement 
when the particle number is large, and hence is less than 
optimal. 



B. Fermionic fields 

Identifying a characteristic COM uncertainty for a cold 
Fermi gas is a natural extension to the above discussion. 
We therefore consider a spin-polarized gas of N identical 
fermions at zero temperature. In other words, wc con- 
sider a system in which all modes below the characteristic 
Fermi energy contain a single fermion each, whilst those 
above this energy are vacant. Using the Hartree-Fock as- 
sumption of wavefunction factorization, this state can be 
written 



ivi/W) = I n 1 10) 



(33) 



where the fermion creation operator for the j-th free- 
particle energy level is defined as 



2] = J d^xx, (x)^i't(x) 



(34) 



and Xj (^) is an orthonormal set of mode functions (en- 
ergy eigenstate basis). 

We now calculate the expectation value of the mean 
position variance operator given by Eq. [551 For the first 
term in this equation, we employ the commutation rela- 



tion 



'I'(x), dj + = Xj(x) to show that 



*(x) 



3=0 



(35) 



where 



'?) 



nt„' (4) 



|0). Taking the 



inner product of this result with it's conjugate leads to 

^' = 1^E / rf''x{|Axnx,(x)|}. (36) 

To evaluate the second term in the expectation value of 
Eq.[2n]we consider the action of two Fermi field operators 
on the system state vector: 



*(y)#(x) 



N 



^(„l).+^-'^(fc-.)(l„J^.,) 
jfe=l 



xxj(x)xfe(y) 



(37) 



Again taking the norm of the result, we find: 

N 

Cixp:) = E d^^d^y{^-y[\xj{^)?\xk{y) 

3k=l 

-Xj(x)xfc(x)xfe(y)xi(y)]} 

Up until now, the single particle basis Xj(^) has been 
arbitrary. For dcfiniteness, we will assume that the func- 
tions Xji'^) — are energy eigenstates of a simple 
harmonic oscillator of mass m and frequency uj. That is. 



(38) 



1 E / rf''x{|xp(j-|x)(x|j)} 

N-l 



3=0 



where we have also defined the first quantized position 
operator q — {q^, %, cjz) such that 9^|x) = Xfj_\x). Recall- 
ing that this operator can also be expressed in terms of 
raising ((5^ ) and lowering (d^) operators, we find 



1 / fi \ ^ 



N \2muj 



and hence the variance due to the wavepacket extent is 
in this case 



Dh 
2mL0 



-N. 



(39) 



The normally ordered component of the COM variance 
can be evaluated using a similar approach. The first line 
in Eq . 1381 vanishes due to the odd parity of the integrand. 
The remaining term gives 

Af-l D 

MXp:) = -EE(-?1'?M|fc)W'?Mli) 

jfe=0 /J=l 



Dh 

2muj 



jk=0 
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which after combining the two terms in brackets leaves 



Dh 
2muj 



N{N - 1). 



(40) 



Taking into account spin degeneracy of an S'-component 
Fermi gas, which gives rise to S independent measure- 
ments each with N/S particles, the COM variance for 
an A^-particle Fermi gas at zero temperature held in an 
harmonic trap is 



lAxP) 



iV2 ■ 



(41) 



The contribution from the normally ordered term van- 
ishes for A*" = 1 and is less than zero for all TV > 1 . As this 
term always vanishes for a bosonic number state, while 
the wavepacket variance is identical for a bosonic gas con- 
strained to the same density profile as the fermionic gas 
we have just considered (i.e. {n{x)) = J2f=i IXj(^)P): 
it is clear that the standard quantum limit for COM 
measurements is intrinsically lower for fermions than for 
bosons. In the case of a harmonic trap, the variance is 
N times smaller for fermions than for bosons. 



IV. QUANTUM WAVE PACKET SPREADING 

In pure state evolution, any increase in the uncertainty 
above the standard quantum limit is regarded as originat- 
ing in quantum wave-packet spreading of the center-of- 
mass. We thus define 



aQM = J (AX2) - ctIq^ 



(42) 



to measure this phenomenon. We note that states for 
which aQM 2> are the "strange quantum states" which 
are mentioned in the photonic case by some previous 
workers;^], and also play a role in the theory of quantum 
non-demolition position measurements and gravity-wave 
detection. This must be carefully distinguished from un- 
certainties in mixed states, which are not caused by quan- 
tum superpositions. 

A similar property to quantum wave-packet spreading 
is already known, and observed experimentally. In the 
Gordon-Haus effect'^ in amplifying optical fibers, the 
soliton time-of-arrival is perturbed by the effect of laser 
amplification. Soliton timing jitter of this type is an ex- 
trinsic property of the amplifying system. Instead, we 
will treat the intrinsic quantum effects which are present 
even in the absence of the spontaneous emission noise of 
a laser amplifier. 

To solve for the time evolution of the COM uncertainty 
of arbitrary states, we first note that Eq.|4]can be directly 
integrated to give the following result, for any (effectively 
non-relativistic) system prepared in an eigenstate of the 
total particle number operator: 

■D 



where m can be either the mass of an atom or the effective 
mass of a polariton, as discussed in section III Al above. 
For extreme relativistic or massless photon propagation, 
the lack of a longitudinal dispersion mechanism implies 
that wave-packet spreading is analogous to diffraction, 
and occurs transversely to the propagation direction. 
This must be treated using the conserved number density 
current operator of Eq. (jlOp . 

This result can be generalized to other states, provided 
that the total number is conserved during the evolution. 
Under this condition, the Heisenberg picture evolution of 
the intensive and quasi-intensive mean position operators 
are respectively 



mN 



and 



Xit) = X(0) +t- 



(43) 



(44) 



In the pulse frame, taking the squares gives the evolution 
of the (true intensive and quasi intensive) mean position: 

(|Ax(t)|2) = (|A(0)|2) + -((i(0).P + P.i(0))/iV) 

m 



+ ^(|P|ViV2), and 



(45) 



(|AX(t)n = (|X(0)n + -(X(0)-P + P-X(0))/(iV) 

m 



+-.{\P\')/{N)' 



Bosonic coherent state evolution 



(46) 



When the system is initially prepared in a coherent 
state (projected to remove the vacuum state, as in Eq. lSOp 
the mean position uncertainties evolve as 



(|AX(i)| = 



^ht fht 

a + 13 h7 — 

m \m 



(47) 



mN 



where a, /9, 7 are functions of the initial wave-packet 
shape only. These are given by: 

a = j |xnx(x)|2d^x 

/3 = ^ y x.[(Vx*(x))x(x)-x*(x)(Vx(x))]d^x 

7 = 1 |Vx(x)pd^x 

where x(x;) is again the normalized spatial mode in the 
co-moving frame. 

This result agrees with previous approximate lin- 
earized results obtained for initial coherent solitonsi4il. 
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but is much more general. It is valid for any initial co- 
herent state, at large photon number, independent of the 
nonlinear interaction, due to the fact that the coefficients 
depend only on the initial coherent state. If we recall 
that the Hamiltonian corresponds to coupled propaga- 
tion of massive bosons, then the reason for this exact 
independence of the coupling is transparent. The defini- 
tion of X corresponds to a center-of-mass measurement, 
which never depends on the two-body coupling. How- 
ever, the standard quantum limit does depend on the 
nonlinearity. This is because the coupling changes the 
pulse-shape (^''^ (x)^'(x)), as an initially coherent pulse 
with uncorrelated bosons evolves into a correlated state. 
Thus, for example, if a ID coherent source produces 
a sech input pulse, so that x{^) — sech(x/2)/2, then 
a = 7r^/3 m^]l3 = 0; 7 = 1/12 m~^. This is not a mini- 
mum uncertainty state, as: 



(AX(0) )p{AP^)c = 7rV36 0.27415.. > 0.25. (48) 



Coherent Gaussian input pulses, on the other hand, are a 
minimum uncertainty state in momentum and position, 
as one might expect. However, these have a large con- 
tinuum radiation when used to form solitons. This leads 
to a paradox: the sech-type coherent soliton has a larger 
uncertainty product than a coherent Gaussian pulse, yet 
experimentally appears more localized! 

More complex input states than coherent states can 
be considered, and they also can be treated exactly. For 
example, Haus and Lai[3| have considered an input state 
of spatially correlated photons, consisting of a superposi- 
tion of eigenstates of the interacting Hamiltonian. In this 
case, the mean pulse-shape is a sech pulse. However, due 
to boson-boson correlations, the resulting state is a mini- 
mum uncertainty state in position and momentum[27] . It 
is unlikely that any existing laser source can produce the 
required correlations. Nevertheless, this example shows 
that, when dealing with correlated states, it is possi- 
ble to reach the minimum uncertainty limit with a non- 
Gaussian pulse-shape. 



V. SOLITON PROPAGATION 

We next consider weakly interacting particle distribu- 
tions confined to a single transverse mode of a waveg- 
uide, which can therefore be treated by an effective one 
dimensional field theory. Examples of such systems in- 
clude photonic wave packets in single-mode optical fibers 
and atomic BEG clouds moving in ID wave guides. Such 
physical systems have the property that they can form 
solitons, in which the classical wave-packet spreading is 
minimized, thus forming excellent candidates for obser- 
vation of these center-of-mass uncertainties. 



A. Effective field tlieory 

Our starting point is the standard Hamiltonian oper- 
ator describing a ID ensemble of spinless bosons which 
interact through a simple delta- function potential. In 
terms of a dimensionless particle density amplitude, (p, 
the Hamiltonian can be written: 



H 



1 



(49) 



where the field operators have the commutation relation 

[0c(r),4(r)] = i<5(C-C')- (50) 

Heisenberg's equations of motion for the field operators 
are then 



(51) 



These equations simultaneously describe at least two 
physically distinct many-body systems, the details of 
which are implicit in the characteristic time to = t/r and 
length a:o = ^/C which link the physical and dimension- 
less coordinate systems. Together with the mean particle 
number n and the sign of the nonlinearity, these scaling 
factors completely specify a particular system. 

Loss (or gain) mechanisms can be formally incorpo- 
rated into this model by including terms in the Hamil- 
tonian which couple the system to external degrees of 
freedom. After tracing over these reservoir states, this 
procedure leads to a master equation for the reduced 
density operator in generalized Lindblad form, which in- 
cludes c-number damping coefficients. Although these 
are frequently added in an ad hoc fashion, it is impor- 
tant to note that in doing so one implicitly assumes that 
the reservoir is at very low temperature and has no ther- 
mal contribution to the system modes. 31 1 



B. Numerical results 

Typical results for the mean position variances dur- 
ing the propagation of a 10'* particle coherent state soli- 
ton with initial shape {(j>(^) ~ ^sech{^() are shown in 
Fig. mja). Here the dashed line represents the analytic 
solution for the total variance (AX^)/x§, while the solid 
lines are numerical results generated using a stochastic 
computer simulation equivalent to the quantum nonlin- 
ear Schrdinger equationfil]. These results include the to- 
tal variance (in complete agreement with the analytical 
result), the standard quantum limit (t^q^ and the mea- 
sure of quantum wave packet spreading Uqj^j. For com- 
parison, Figs.[2l^b) and[2l[c) display the evolution of these 
variances for wider and narrower pulses than the classical 
soliton solution. 

The standard quantum limit for the soliton is constant, 
which indicates each soliton envelope remains nearly in- 
variant; as one might expect from a classical analysis. 
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Figure 2: Mean position variances for ID coherent pulses of 
n = 10* particles propagating in a dispersive medium. The 
three figures correspond to three different wavepacket shapes: 
(a) the soliton solution to the ID nonlinear Schrodinger equa- 
tion (0(C)) = |sech(iC)i (b) a wider pulse of the form 
{<i>{C)) oc sech(iC) and (c) a narrower pulse (<^(C)) oc sech(C). 
The solid lines represents the analytically determined total 
variance {AX'^) /xq, while the dashed and dash-dotted lines 
represent the numerical results for the classical (csqi^/xq) 
and quantum {(^qm/xo) contributions respectively. 



This is a remarkable macroscopic quantum effect. An 
optical coherent state soliton consists of linear super- 
positions of continuum photons and different number- 
momentum eigenstates. This is clearly different to a 
classical soliton with frequency jitter, yet it behaves very 
similarly as far as the center position is concerned. On 
the other hand, the standard quantum limit for a pulse in 
a linear medium increases in the same way as the soliton 
quantum wave-packet spreading, due to the dispersive 
spreading in the average intensity. In the linear case, an 
initial coherent state remains coherent, so the increased 
quantum positional uncertainty can be attributed to the 
shot-noise error intrinsic to the detection of a pulse whose 
envelope is not sharply localized. 

This indicates that the intrinsic quantum wave-packet 
spreading can be observed only for the soliton, where it 
is easily distinguishable from coherent shot-noise effects. 
However, this does not mean that a quantum soliton is 
unstable. Each soliton preserves its soliton pulse shape. 
In this sense, the soliton quantum diffusion is different 
from the classical pulse broadening due to linear disper- 
sion. In fact, in the case of the correlated initial state, 
the second-order and the fourth-order correlation func- 
tions are invariant 1291. 



VI. PRACTICAL EXAMPLES 

In this section we consider some practical examples 
and numerical estimates with currently available experi- 
mental technologies. 

A. Photonic systems 

In the optical case, the soliton propagation model de- 
scribes a distribution of interacting photons (or strictly, 
polaritons) propagating in a single-mode fiber with a 
third-order Kerr nonlinearity, x*''^^- Here the sign of the 
interaction term is positive for normal dispersion and 
negative for anomalous dispersion. For x^"^) soliton prop- 
agation, we require the carrier frequency to be inside the 
anomalous dispersion regime. 

Typical orders of magnitude of soliton parameters are 
n = 10^, to = 10~^^s, xq = Im. With these parameters, 
we obtain a standard quantum limit of 10~*m for a single 
coherent pulse. In 1km of propagation, the wave packet 
would therefore spread to 10~^m, much greater than the 
standard quantum limit. Although still less than the ac- 
tual soliton extent of lO^^m, this is a remarkably large 
quantum wave-packet spreading in a composite object 
of 10^ particles. An experiment on quantum-mechanical 
wave-packet spreading in a composite object of this type 
(such as that described in would test quantum me- 
chanics in a region of much larger particle number than 
previously achieved. 

Unfortunately, however, even fiber attenuations lower 
than the usual O.ldB/km translate to losses of over 10^ 
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photons from such a pulse over this distance. Since each 
particle lost can be used to obtain information about 
the mean position of a pulse, these losses would destroy 
the coherence properties of a wave packet. Even so, the 
remaining statistical uncertainty is a practical concern - 
it is, for example, likely to have a severe effect on the 
error-rate in pulse-position logic [28|. The effect may be 
enhanced by using dispersion-engineered fibers 



B. Degenerate Bose gasses 

The same effective field theory can also describe a di- 
lute, single-species gas of massive bosonic particles, in- 
teracting through low-energy S-wave (/ = 0) scattering 
events and held in a potential which confines the parti- 
cles to a single transverse mode but allows longitudinal 
propagation. In this case, one finds that: 



to 



mngiD 



(52) 



where m is the single-particle mass, andoiD = 2hasUj± 
is the effective ID interaction strength [30|. For temper- 
atures well below the BEC transition temperature, such 
gasses are usually considered to be the matter- wave ana- 
log of coherent laser output, and are thus sometimes re- 
ferred to as atom-lasers. 

Although there are various paths one can take to real- 
izing matter- wave solitons, we consider here only bright 
solitons in the absence of any periodic potential - i.e., 
solitons for which the required nonlinearity is provided 
by an attractive interaction between atoms. Such solitons 
have been observed [13, [1^ in BECs of ''Li atoms, where 
the presence of a Feshbach resonance is used to tune the 
effective scattering length from repulsive to attractive. 
These experiments present the very real possibility of di- 
rectly observing the quantum dispersion of a mesoscopic 
object composed of massive particles, thanks to the ex- 
tremely slow rate of atom loss from these systems . For 
instance, the total predicted loss rate for a trapped cloud 
of attractive ^Li atoms of roughly the same number as 
that in [s^ is around 50/s.[33| (This is analogous to an 
optical pulse traveling along a fiber with an attenuation 
of less than IQ-MB/km!) 

Drawbacks in dealing with atomic clouds include the 
effect of finite temperatures. This places a classical un- 
certainty in the center of mass momentum and therefore 
leads to a classical spread in center of mass position over 
time, obscuring the quantum diffusion. For a ID gas of N 
bosons of mass m and temperature T, the equipartition 
theorem gives 



crth 



mN 



(53) 



where &th is the linear rate of increase in the thermal 
mean position uncertainty. 



For example, using numbers from the experiment of 
Khaykovich et al. [32l | and assuming a condensate temper- 
ature of 0.5rc ~ 80nK we find that the standard devia- 
tion in the mean position due to thermal noise increases 
linearly at a rate of 1.3 x 10~''m/s. This is quite fast, 
when compared to the effect of quantum spreading of the 
soliton wave packet, which increases the standard devia- 
tion at the much slower rate of 4.2 x lO^^m/s. As this is 
an order of magnitude smaller than the thermal uncer- 
tainty, whose growth rate goes as the square root of the 
cloud temperature, this implies that temperatures less 
than InK are required to directly observe quantum wave 
packet spreading. Although this is well below what is 
achieved in current ^Li experiments, Leanhardt et al. [35| 
have achieved temperatures of less than 0.5nK in con- 
densates of spin-polarized ^'^Na, indicating that it may 
not be impossible to achieve such low temperatures in 
matter-wave soliton experiments. We note that adiabatic 
changes in the interaction strength could be used to com- 
press the soliton, thereby increasing the quantum posi- 
tion uncertainty relative to the standard quantum limit. 



C. Degenerate Fermi gasses 

Ultra-cold degenerate Fermi gase s have recently been 
obtained in several laboratories]!^ [s^]. The prospects 
for fermionic solitons have not yet been established 
clearly, although tunable Feshbach resonances with both 
attractive and repulsive interactions are known to exist. 
The species observed experimentally include ""^if , ^Li 
and ^He*. Atom correlations ^SSj have already been ex- 
perimentally observed in optical measurements with ^'^K. 
One of the most interesting cases is metastable fermionic 
Helium, in which atoms can be counted directly using 
multi-channel plate detectors [131 . This affords an in- 
teresting possible avenue to testing the prediction that 
center-of-mass fluctuations are reduced for fermions rel- 
ative to bosons. 



D. Double-slit interference of soliton 'particles' 

Provided one could overcome this difhculty, it would be 
extremely interesting to conduct double slit interference 
experiments with such large quantum objects. One pos- 
sible implementation of this idea - essentially an atomic 
implementation of the fiber experiment proposal in [^ - 
is shown in schematic as Fig. [3] Here a matter- wave soli- 
ton is allowed to propagate without disturbance until its 
quantum wave packet is large compared with the soliton 
width. A tightly focused laser pulse is then used to re- 
move the central and outer elements of the distribution, 
effectively creating a pair of apertures through which the 
remaining components of the distribution pass. At some 
later time, the distribution is imaged and the mean po- 
sition of the soliton measured. Another approach - pro- 
vided the separation was large compared to the soliton 
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Figure 3: Possible soliton double slit interference experiment. 
The soliton wave packet (I) broadens over time due to the 
linear increase in the quantum uncertainty in the mean posi- 
tion. Lasers then eliminate solitons outside of two 'slits' (II). 
Any remaining soliton is left to propagate until it undergoes 
absorption imaging (III). 

size - would be to use temporal switching of a localized 
reflective potential, in order to only affect one part of 
the quantum soliton superposition. Averaging over many 
shots should then produce an interference pattern with 
any thermal effects being manifest as a loss in fringe vis- 
ibility. 

One of the primary challenges one would face in con- 
ducting this experiment would be addressing the fact that 
its duration needs to be less than the average particle 
lifetime in order to maintain coherence, while still long 
enough to allow sufficient quantum wave packet spread- 
ing. A possible solution might be to use a blue-detuned 
laser to provide a small repulsive potential near the cen- 
ter of the wave packet to enhance the broadening and 
thereby reduce the length of time required to produce 
an interference pattern. This technique might also al- 
low one to forgo the use of destructive laser pulses to 
create the apertures by instead relying on this repulsive 



potential to separate the center of mass wave packet into 
two components. Another technique which appears pos- 
sible, is to employ the recently proposed soliton quantum 
beam-splitter |39| to initially separate two quantum soli- 
ton components based on their velocities, prior to subse- 
quent recombination. 

A related challenge would be resolving the resulting 
fringes, as their width is inversely proportional to parti- 
cle number and limited by the trap lifetime through the 
maximum length of the experiment. (This can be shown 
by treating the soliton as a de Broglie wave of mass miV; 
m and N being the single-particle mass and the total 
number of constituent atoms respectively.) 

Despite these issues, such an interference pattern 
would provide smoking gun evidence that quantum me- 
chanical superpositions of massive composite objects of 
mesoscopic scale had been achieved. 

VII. SUMMARY 

In summary, quantum wave-packet spreading is a re- 
markable macroscopic quantum effect. It provides a fun- 
damental limit - independent of amplifier noise - both 
to high-speed communications in dispersive waveguides, 
and to applications of atom lasers. Since it provides a 
means by which quantum mechanics can couple to the 
gravitational field, it may also provide a route to new 
tests of quantum mechanics. Of more practical interest 
is the fact that the standard quantum limit for center-of- 
mass measurement variance is highly sensitive to particle 
statistics. It is decreased by a factor of TV when the par- 
ticle statistics are fermionic. This appears testable in 
atom-counting experiments with meta-stable ^He . 
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